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We present the complete calculation of W-boson production in association with a jet in hadronic 
collisions through next-to-next-to-leading order in perturbative QCD. To cancel infrared divergences 
we discuss a new subtraction method that exploits the fact that the N-jettiness event-shape variable 
fully captures the singularity structure of QCD amplitudes with final-state partons. This method 
holds for processes with an arbitrary number of jets, and is easily implemented into existing frame¬ 
works for higher-order calculations. We present initial phenomenological results for W+jet produc¬ 
tion at the LHC. The NNLO corrections are small and lead to a significantly reduced theoretical 
error, opening the door to precision measurements in the W+jet channel at the LHC. 


INTRODUCTION 

Run I of the Large Hadron Collider (LHC) was a re¬ 
markable success, culminating in the 2012 discovery of 
the long-awaited Higgs boson. One major contributor to 
this exciting outcome was the precision QCD framework 
used to model and understand both the sought-after sig¬ 
nals and the often overwhelming backgrounds. The de¬ 
mand for ever-more sophisticated precision QCD calcu¬ 
lations is only increasing as Run II of the LHC begins. 

Calculations through next-to-next-to-leading order 
(NNLO) in perturbative QCD are becoming increasingly 
necessary to match the precision of LHC measurements. 
During the past years the first few NNLO calculations 
for 2 -+ 2 scattering processes, needed for precision 
phenomenology in many channels at the LHC, have be¬ 
gun to appear. The full results for tt production [1], 
single-top production [2] and for several color-neutral fi¬ 
nal states are available [3] . The situation is less advanced 
for processes containing final-state jets, which possess a 
more complex singularity structure. Partial results are 
available for inclusive jet production [4] and Higgs+jet 
production [5]. All such NNLO calculations require a 
subtraction scheme to extract infrared singularities from 
real-emission amplitudes and cancel them against the 
corresponding divergences appearing in the virtual am¬ 
plitudes. Currently, two subtraction schemes have been 
demonstrated to successfully handle processes contain¬ 
ing final-state jets in hadronic collisions: antennae sub¬ 
traction [6] and sector-improved residue subtraction [7]. 
Although powerful, both approaches require the devel¬ 
opment of a sophisticated and specialized machinery. It 
appears difficult to merge such calculations into a com¬ 
mon framework with other standard tools used by the 
high energy physics community. 

An alternative to such techniques is ^T-subtraction [8] . 
It leverages the observation that for color-neutral final 
states, the transverse momentum qj- of the final state 


completely determines the singularity structure of the 
contributing QCD amplitudes. For qp above some small 
resolution parameter, the result is simply the NLO cal¬ 
culation for the color-neutral state plus an additional jet. 
Below the cutoff, the cross section is dominated by large 
logarithms of the final-state mass scale over qp. This 
cross section can be obtained by expansion of the well- 
known CSS resummation formula [9] to fixed order. Un¬ 
fortunately, ^-subtraction cannot be extended beyond 
color-neutral final states, since the transverse momentum 
no longer completely describes the singularity structure 
of QCD cross sections with final-state jets. 

In this manuscript we discuss a subtraction scheme 
that overcomes these various limitations. It is valid for an 
arbitrary number of final-state jets, and maximally reuses 
information available in existing NLO calculations. It is 
based on the observation that the X-jettiness event shape 
variable Tn [10] controls the singularity structure of QCD 
cross sections with N final-state jets. One can partition 
the phase space according to a resolution parameter Tfj ut . 
Above the cutoff, the result is the NLO calculation with 
N + 1 jets. Below the cutoff, the resummation formula 
for Tn can be used to obtain the cross section through 
NNLO. In this phase-space region the cross section can 
be written in terms of a small number of known universal 
functions, together with the process-dependent two-loop 
virtual corrections. A similar division of the final-state 
phase space using the invariant mass of the hadronic radi¬ 
ation was used to obtain top-quark decay at NNLO [11]; 
we discuss here the general applicability of this idea to 
X-jet cross sections through the use of Tn- 

We demonstrate our method with a highly non-trivial 
example: W+jet at NNLO. This calculation is needed for 
numerous LHC phenomenological applications, including 
the determination of the gluon distribution function. We 
present here selected numerical results for W+jet produc¬ 
tion for the LHC. We find that the NNLO corrections de¬ 
crease the NLO cross section by approximately 1%, and 
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significantly reduce the residual theoretical uncertainty 
as estimated by scale variation. Our computation makes 
precision measurements possible in the W+ jet process at 
the LHC. 


THEORETICAL FRAMEWORK 

We sketch here the construction of the A-jettiness sub¬ 
traction scheme. We begin with the definition of N- 
jettiness,7/v, a global event shape variable designed to 
veto final-state jets [10]: 

The subscript N denotes the number of jets desired in the 
final state, and is an input to the measurement. For the 
W+jet process considered here, we have N = 1. Values 
of 71 near zero indicate a final state containing a single 
narrow energy deposition, while larger values denote a 
final state containing two or more well-separated energy 
depositions. The pi are light-like vectors for each of the 
initial beams and final-state jets in the problem, while the 
qk denote the four-momenta of any final-state radiation. 
The Qi are dimensionful variables that characterize the 
hardness of the beam-jets and final-state jets. We set 
Qi = 2 Ei, twice the energy of each jet. 

To proceed, we first note that at NNLO, the cross sec¬ 
tion consists of contributions with Born-level kinematics, 
and processes with either one or two additional partons 
radiated. We partition the phase space for each of these 
terms into regions above and below T§ ut '. 

&nnlo= J d$ JV |M JV | 2 + J ddw+i \Mn+i\ 2 6n 

+ J d<fw+2 \Mn+2\ 2 6n + J d'hjv+i |.Adjv-i-i | 2 9 ^ 

+ J d<lw +2 \Mn+2\ 2 

= &NNLO (Tn < T^ Ut ) + &NNLo(Tn > T^) 

( 2 ) 

Here, we have abbreviated 9 ^ = 9(T£ ut — Tn) and 
Tv = 9(Tn — T£ ut ), and have suppressed for simplicity 
the allowed introduction of any infrared-safe measure¬ 
ment function under the phase-space integral. The first 
three terms in this expression all have T/v < T^ ut , and 
have been collectively denoted as ctnnlo(Tn < T£ ut )- 
The remaining two terms have Tn > Tjf , and have been 
collectively denoted as <Jnnlo(Tn > T£ ut )- Contribu¬ 
tions with Born-level kinematics necessarily have Tn = 0. 
We note that a similar partitioning of phase space to 
separate the singular and non-singular regions in Tn has 
proven quite useful in the context of merging fixed-order 
calculations with parton showers in the effective-theory 
framework pioneered by the Geneva collaboration [12]. 


The key advance that allows us to compute the cross 
section to NNLO below 7^ ut is the existence of a fac¬ 
torization theorem that gives an all-orders description of 
fV-jettiness for small Tn [13, 14], Using this factorization 
theorem, the cross section for a hadronic process with Tn 
less than some value Tfi ut can be written in the schematic 
form 


cr(T N <TN Ut ) = J 


N 


IP- 


(3) 


Here, H describes the effect of hard radiation; when di¬ 
mensional regularization is used this function simply en¬ 
codes the virtual corrections to the process. B encodes 
the effect of radiation collinear to one of the two ini¬ 
tial beam directions; the importance of the beam func¬ 
tion more generally in describing hadronic collisions was 
first pointed out in Ref. [14]. It can be further decom¬ 
posed as a perturbative matching coefficient convolved 
with the usual parton distribution function. S describes 
the soft radiation, and J n contains the radiation collinear 
to a final-state jet. The ellipsis denotes power-suppressed 
terms which become negligible for Tn •C Qi ■ The deriva¬ 
tion of this all-orders expression in the small-7/v limit 
relies heavily upon the machinery of Soft-Collinear Ef¬ 
fective Theory [15]. 

If this formula is expanded to fixed-order in the strong 
coupling constant, it reproduces the fixed-order cross sec¬ 
tion (Jnnlo(Tn < T^ ut ) for low 7needed in Eq. (2). 
Currently the hard function is known at NNLO for 
many phenomenologically interesting cases, including the 
W+jet process considered here [16]. The beam functions 
are known at NNLO [17], as are the jet functions [18] and 
the soft function [19]. 

A full NNLO calculation requires as well the high Tn 
region above T^ ut . However, a finite value of Tn implies 
that there are actually N + 1 resolved partons in the 
final state. This is the crucial observation; Tn completely 
describes the singularity structure of QCD amplitudes 
that contain N final-state partons at leading order. The 
high Tn region of phase space is therefore described by a 
NLO calculation with N + 1 jets. We must choose T£j ut 
much smaller than any other kinematical invariant in the 
problem in order to avoid power corrections to Eq. (3) 
below the cutoff. We address this issue explicitly in the 
context of W+jet in the next section. 

To summarize, we list here the exact steps needed for 
our calculation of the full NNLO calculation of W+jet. 


• Generate an event for the N + 1-jet process at 
NLO, which may contain either N + 1 partons 
or N + 2 partons. Determine the reference vec¬ 
tors pi in Eq. (1) by performing a pre-clustering of 
the radiation using a jet algorithm, as discussed in 
Refs. [10, 20]. The determination of the Pi is insen¬ 
sitive to the choice of jet algorithm in the small- 
T^ ut limit [10]. 
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• Calculate Tn according to Eq. (1). If Tn > T^ ut , 
keep the event. Events satisfying this criterion form 
the NLO cross section for the N -f 1-jet process 
&nnlo(Tn > T^ ut ) needed in Eq. (2). This NLO 
cross section can be obtained using any standard 
technique or code. If Tn < T^ ut , reject the event. 

• Obtain the cross section ct(Tn < T§ ut ) by expand¬ 
ing the resummation formulae of Eq. (3) to NNLO. 
We note that this fully accounts for all V-parton 
contributions in Eq. (2), as well as for terms with 
additional partons where Tn < T^ ut . 


VALIDATION OF THE FORMALISM 



TT [GeV] 


We next discuss the validation of the T/v-subtraction 
formalism, in the context of our calculation of W+jet at 
NNLO. An advantage of our formalism is that it max¬ 
imally reuses known information coming from existing 
NLO calculations. Above T^ ut we need a NLO calcula¬ 
tion of TV+2-jets, which we obtain from MCFMv8.0 [21]. 
For the terms which contribute below Tfi ut , we have 
checked our implementation of the two-loop virtual cor¬ 
rections against those contained in PeTeR [22]. Our cal¬ 
culation and validation of the necessary V-jettiness soft 
function has been detailed in a separate publication [19]. 

A powerful check of our formalism is that the loga¬ 
rithmic dependence on Tfi ut that occurs in the separate 
low and high Tn regions cancels when they are summed. 
We show in Fig. 1 the results of this check. The plot 
shows only the 0(a correction to the cross section as a 
function of T^ ut ', we have checked that the O(a^) NLO 
cross section obtained with this technique agrees exactly 
with the known results. These cross sections are obtained 
using CT14 parton distribution functions [23], and con¬ 
tain the following fiducial cuts on the final-state jet from 
CMS [25]: p> 30 GeV, \rj jet \ < 2.4. The ATLAS anal¬ 
ysis is similar but with slightly different cuts [26]. Both 
the renormalization and factorization scales have been 
set to p = Mw and varied from this choice by a factor of 
two. Over the region 0.05 GeV < T^ ut < 0.08 GeV, the 
sum is stable to better than one picobarn, a size which 
represents an 0.1% correction to the total cross section. 
The NLO corrections to W+2-jets, have been obtained 
using the double precision version of MCFM, both the 
single-core version and the new multi-core implementa¬ 
tion [24]. We have checked that for larger values of T^ ut , 
the power corrections in Eq. (3) begin to become impor¬ 
tant. For the numerical results presented in the remain¬ 
der of this paper we use the choices T£ ut = 0.05, 0.06 
GeV to cross-check their Tfi ut independence. 

As a final check of our computation, we have applied 
our formalism to calculate the NNLO corrections to Higgs 
production in association with a jet, for which partial 
results are available [5]. We find agreement with the 


Figure 1: The NNLO correction Aonnlo normalized to the 
NLO cross section as a function of Tf ut ■ The solid blue line 
denotes the result for the central scale choice p = Mw, while 
the red hatched area indicates the correction in the range 
Mw / 2 < p < 2 Mw- The solid cyan and green lines indicate 
the results for p = 2Mw and p = Mw/2, respectively. The 
NNLO result for each scale choice has been normalized to the 
NLO cross section obtained using the same scale choice. 

fiducial cross section for this process. The details of this 
computation, together with phenomenological results for 
the LHC, will be presented in a separate manuscript [27]. 

NUMERICAL RESULTS 

We present here first numerical results for W + +jet 
production at the LHC. We focus on y/s = 8 TeV col¬ 
lisions, and use CT14 parton distribution functions [23] 
at the same order of perturbation theory as the corre¬ 
sponding partonic cross section. Our central scale choice 
is p = M\v- To obtain an estimate of the theoretical er¬ 
rors we vary p away from this choice by a factor of two. 
We again impose the following cuts on the final-state jet: 

> 30 GeV, \'f]j e t\ < 2.4. We reconstruct jets using 
the anti-fey algorithm [28] with R = 0.5. 

We begin by showing in Table I the total cross sec¬ 
tion subject to the above cuts at LO, NLO, and NNLO 
in the strong coupling constant. These numbers include 
the branching fraction for the W-boson to decay to a 
single lepton flavor. We note that the numerical error 
on these numbers is at the several-per-mille level. The 
cross section shifts by +40% when going from LO to NLO 
in perturbation theory, but only by approximately +1% 
when going from NLO to NNLO. The scale variation is 
approximately ±7% at LO and NLO, while at NNLO it 
is reduced to the percent level. We note that at NNLO 
the largest cross section is obtained for p = M w , leading 
to the lack of scale variation in the upper direction in the 
Table below. The residual theoretical error is reduced to 
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the percent level at NNLO, and excellent convergence of 
the perturbative series is obtained. 


p> 30 GeV, | Tjjet | < 2.4 

Leading order: 

56lia5 pb 

Next-to-leading order: 

796l®g pb 

Next-to-next-to-leading order: 

807t? 0 pb 


Table I: Fiducial cross sections, defined by > 30 GeV, 
\r)j e t\ < 2.4, using CT14 PDFs at each order of perturbation 
theory. 

In Fig. 2 we show the transverse momentum spectrum 
of the leading jet at LO, NLO and NNLO in perturba¬ 
tion theory. The ratios of the NLO cross section over 
the LO result, as well as the NNLO cross section over 
the NLO one, are shown in the lower inset. The shaded 
bands in the upper inset indicate the theoretical errors 
at each order estimated by varying the renormalization 
and factorization scales by a factor of two around their 
central value, as do the vertical error bars in the lower 
inset. In the lower inset we have shown the results for 
both T^ ut = 0.05 GeV and T^ ut = 0.06 GeV, for the scale 
choice p = Mw, to demonstrate the T^ ut independence 
in every bin studied. The NLO corrections are large and 
positive for this scale choice, increasing the cross section 
by 40% at = 40 GeV and by nearly a factor of two 
at jpj f* = 180 GeV. The scale variation at NLO reaches 
approximately ±20% for = 180 GeV. The shift when 
going from NLO to NNLO is much more mild, giving only 
a percent-level decrease of the cross section that varies 
only slightly as pi!^ 1 is increased. The scale variation at 
NNLO is at the percent level and is nearly invisible on 
this plot. 

The transverse momentum spectrum of the W-boson is 
shown in Fig. 3. The NLO corrections are again 40% for 
Pr > 50 GeV with a sizable scale dependence, while the 
NNLO corrections are flat in this region and decrease 
the cross section by a small amount. The phase-space 
region p^ < 30 GeV only opens up at NLO, leading 
to a different pattern of corrections for these transverse 
momentum values. The instability of the perturbative 
series in the bins closest to the boundary pjf = 30 GeV 
is caused by the well-known Sudakov-shoulder effect [29] . 


CONCLUSIONS 

We have presented in this manuscript the complete 
NNLO calculation of W-boson production in association 
with a jet in hadronic collisions. To perform this compu¬ 
tation we have discussed a new subtraction scheme based 
on the IV-jettiness event-shape variable Tn ■ We have 
validated our approach in several ways: when possible 
the various components have been cross-checked against 
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Figure 2: The transverse momentum spectrum of the lead¬ 
ing jet at LO, NLO and NNLO in perturbation theory. The 
bands indicate the estimated theoretical error. The lower in¬ 
set shows the ratios of the NLO over the LO cross section, 
and the NNLO over the NLO cross section. The red vertical 
error bars in the lower inset indicate the scale-variation error. 
The blue and black lines in the lower inset respectively show 
the distribution for T^ 1 = 0.05 GeV and T^ ut = 0.06 GeV, 
for the scale choice p = Mw- 

known results in the literature, the necessary cancella¬ 
tion of the logarithmic T^ ut between the phase-space re¬ 
gions Tn > T£ ut and Tn < T^ ut has been established, 
and we have reproduced known results for Higgs pro¬ 
duction in association with a jet at NNLO. The NNLO 
corrections to the W+jet process indicate a remarkably 
stable perturbative series ready to be used for precision 
measurements at the LHC. We will further study the 
phenomenological impact of our NNLO result in future 
work, including the prediction for the exclusive one-jet 
bin, where an intricate interplay between various sources 
of higher-order corrections was recently pointed out [30] . 

We believe that the development of the jettiness- 
subtraction technique represents a watershed moment in 
the field of higher-order calculations. For the first time a 
completely general subtraction scheme valid for any num¬ 
ber of final-state jets has been introduced that is based 
on a factorization theorem which extends to all orders in 
perturbation theory and is straightforward to implement 
in existing frameworks for NLO calculations. We antic¬ 
ipate that the W+jet process presented here is only the 
first of many results obtained with this novel technique. 
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Figure 3: The transverse momentum spectrum of the W- 
boson at LO, NLO and NNLO in perturbation theory. The 
bands indicate the estimated theoretical error. The lower in¬ 
set shows the ratios of the NLO over the LO cross section, 
and the NNLO over the NLO cross section. The red vertical 
error bars in the lower inset indicate the scale-variation error. 
The blue and black lines in the lower inset respectively show 
the distribution for T£ ut = 0.05 GeV and T/v“ 4 = 0.06 GeV, 
for the scale choice fi = Mw- 
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